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ABSTRACT 

The motion of two planets around a Sun-like star under the combined effects of mutual 
interaction and tidal dissipation is investigated. The secular behaviour of the system is ana- 
lyzed using two different approaches. First, we solve the exact equations of motion through 
the numerical simulation of the system evolution. In addition to the orbital decay and circu- 
larization, we show that the final configuration of the system is affected by the shrink of the 
inner orbit. Our second approach consist in the analysis of the stationary solutions of mean 
equations of motion based on a Hamiltonian formalism. 

We consider the case of a hot super-Earth planet with a more massive outer companion. 
As a real example, the CoRoT-7 system is analyzed solving the exact and mean equations 
of motion. The star-planet tidal interaction produces orbital decay and circularization of the 
orbit of CoRoT-7b. In addition, the long-term tidal evolution is such that the eccentricity of 
CoRoT-7c is also circularized and a pair of final circular orbits is obtained. A curve in the 
space of eccentricities can be constructed through the computation of stationary solutions 
of mean equations including dissipation. The application to CoRoT-7 system shows that the 
stationary curve agrees with the result of numerical simulations of exact equations. 

A similar investigation performed in a super-Earth- Jupiter two-planet system shows that 
the doubly circular state is accelerated when there is a significant orbital migration of the inner 
planet, in comparison with previous results were migration is neglected. 

Key words: celestial mechanics - planetary systems. 



1 INTRODUCTION 

It is well known that close-in planets undergo tidal interactions 
with their host stars. The tidal effect produces orbital decay, 
circularization and spin-orbit synchronization of a planet orbit- 
ing a slow-rotating star (Dobbs-Dixon, Lin & Mardling 2004; 
FeiTaz-Mello, Rodriguez & Hussmann 2008; Jackson, Greenberg 
& Barnes 2008). When the tidally affected planet has an eccentric 
companion, its evolution is more complicated. As we shall show in 
this paper for two different systems, new features appear due to the 
interplay of tidal effects and mutual interactions between planets. 

Several works have already investigated the evolution of two- 
planet systems accounting for tidal dissipation (Wu & Goldreich 
2002; Mardling & Lin 2004; Mardling 2007; Zhou & Lin 2008, 
Mardling 2010, Greenberg & Van Laerhoven 201 1). Particularly in 
Mardling (2007), a secular model based on the Legendre expan- 
sion of the disturbing forces up to fourth order in semi-major axes 
ratio, a\/a2, was introduced. She has shown that the outer com- 



panion excites the eccentricity of the inner planet, accelerating its 
migration toward the star due to tidal dissipation. Moreover, the 
evolution of the system under tidal effect follows several stages, 
depending on the circulation or oscillation of the angle between 
the lines of apses. The excitation of the inner planet eccentricity 
reaches a quasi-equilibrium value and, in a long time-scale (which 
depends on physical parameters and initial orbital configuration), 
both eccentricities are damped to zero. If the companion planet can 
sustain a non-zero eccentricity of the inner planet for times of the 
order of the age of the system, the mechanism could explain some 
observed oversized planets, since the body can inflate in response to 
tidal heat. The companion planet hypothesis would thus explain the 
large size of some planets, e.g. HD 209458b, HAT-Plb and WASP- 
Ib (see Mardling 2007 for details). 

Although Mardling's model has no restrictions on the values 
of the planetary mass ratio, its applications have given priority to 
systems with more massive inner planet^. In this case, in a first 
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^ Except in the case of HAT-P-13 non-coplanar system (see Mardling 
2010) 
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approximation, we may assume that the planetaiy semi-major axes 
do not vary during the tidal evolution. This approximation may be 
sufficient for systems with a less massive outer planet, in which 
case the excitation of the inner planet eccentricity is not strong 
enough to induce its rapid orbital decay due to the tidal friction. 
As a consequence, the orbital decay of the inner planet occurs in 
time-scales much longer than the ages of the systems analyzed in 
Mardling (2007). In this scenario, analytical results have shown a 
good agreement with results of numerical simulations, performed 
through numerical integrations of the Lagrange's planetary equa- 
tions. 

In the present work, we use a different approach, with no re- 
strictions on the semi-major axes variation, which allows us to in- 
vestigate the tidal evolution of two-planet systems with an arbi- 
traiy mass ratio. We perform numerical simulations, using New- 
ton's equations of the coplanar motion of two planets and includ- 
ing the forces due to tidal friction according to Mignard (1979). We 
show that the angular momentum conservation constrain the vari- 
ation of the planetary orbital elements, even in the case of rapid 
orbital decay. 

We perform an adaptation of the conservative analytical model 
of secular motion to the systems with slow loss of the energy, such 
as tidally affected pair of planets. Assuming that the planets are far 
enough from any mean-motion resonances, only secular terms are 
maintained in the disturbing function. Previous works have shown 
that for a given value of the total angular momentum, all possi- 
ble motions of the secular system are constant-amplitude oscilla- 
tions around one of the two equilibria of the secular Hamiltonian, 
known as Mode I and Mode II (Michtchenko & Ferraz-Mello 2001, 
Michtchenko & Malhotra 2004). We also know that the location of 
the stationary solutions in the phase space is uniquely defined by 
two parameters: the ratios of planet masses and semi-major axes. 

When we consider the non-conservative secular system, in 
which the rate of dissipation is sufficiently slow (more precisely, 
it is slower than the proper frequency of the system), the parame- 
ter a\/a2 varies slightly due to the slow orbital decay of the inner 
planet. For all possible values of a\/a2 and the total angular mo- 
mentum, we may calculate the locations of the equilibria in the 
phase space. A curve of stationary solutions can be constructed and 
then compared with the result of the numerical simulation of the 
system (see Michtchenko & Rodriguez 201 1 for a detailed discus- 
sion). 

In this paper, we restrict our investigations to systems with a 
hot super-Earth and a more massive outer companion. By super- 
Earth, we mean planets with an upper limit of 10 m®, where 
is the mass of the Earth (Fortney, Marley & Barnes 2007; Valencia, 
Sasselov & O'Connell 2007). Up to date, twenty-four known super- 
Earth have been discovered with masses ranging from 1.94 to 
9.22 mm and the majority of these planets evolve in multi-planet 
systemo Almost half of the known super-Earth have orbital peri- 
ods smaller than 5 days; so small periods indicate that the planets 
have experienced an orbital decay due to the tidal interaction with 
the star. 

We apply our approach to the CoRoT-7 extrasolar system, 
which is composed by one super-Earth and a more massive outer 
companion (Queloz et al. 2009; Leger at al. 2009; Ferraz-Mello 
et al. 2011). Previous investigations involving tidal evolution of 
CoRoT-7 planets include the analysis of tidal heating of CoRoT- 
7b and the consequences for its orbital migration. It was shown 
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that the past history of the planet may have included a tidal heating 
large enough to be compared with the corresponding tidal heating 
of the Jovian satellite lo (Barnes et al. 2010). The coupled interac- 
tion between tides and secular evolution in the CoRoT-7 system is 
investigated in this paper showing that a final orbital configuration 
with two circular orbits is obtained in a few Myr. 

This paper is organized as follow: In Section|2]we briefly de- 
scribe the tidal evolution for a single-planet system. In Section [3] 
we analyze the two-planet secular behaviour through the numer- 
ical integrations of the exact equations of motion including the 
tidal force (in the coplanar case). We consider two applications: 
the real system CoRoT-7 and an hypothetical system containing a 
hot super-Earth with a Jupiter companion. The angular momentum 
conservation of the system is discussed in Section |4] providing a 
simple explanation for the damping of the outer planet eccentricity 
and the orbital decay of the inner planet. In Section |5] we inves- 
tigate the secular dynamics of CoRoT-7 planets through the study 
of stationary solutions of the mean equations using the Hamilto- 
nian formalism. The role of stellar tides is analyzed in Section [6] 
Finally, Section [7] is devoted to provide a general discussion and 
conclusions. 



2 TIDAL EVOLUTION OF SINGLE-PLANET 

We first consider the tidal interaction between one planet and its 
host star. We call mo, Ro, the mass and stellar radius, whereas mi, 
R\ are the cortesponding planetary values. In the case of zero or- 
bital inclination, the mean variations of orbital elements produced 
by the joint effect of both planetary and stellar tides are given by 



<h>= --na~^s[{l + 23e^) + 7e^D] 
and 



(1) 



<e>^--nea~^s[9 + 7D], (2) 

where D = p/2s and 

9 fco mi 5 , 9 ki mo 5 
4 Qo JTiQ 2 Qi mi 

are the strengths of the stellar and planetary tides, respectively 
(Dobbs-Dixon et al. 2004; Ferraz-Mello et al. 2008). In above equa- 
tions, a, e and n are semi-major axis, eccentricity and mean orbital 
motion, respectively, ki and Qi are the Love number and dissi- 
pation function of the deformed body, for i = 0, 1. The dissipa- 
tion function is associated to the efficiency for which the energy 
is dissipated in the interior of the tidally deformed body. The re- 
sults are valid up to second order in eccentricity, assuming a quasi- 
synchronous planetary rotation and a slow-rotating star (see Section 
I3.2.3I I. To obtain the variations produced by stellar tide, a linear 
model with constant time lag was assumed to relate phase lags and 
tidal frequencies (see Section [3Tt . 

From equations ([T) and (|2j we obtain da/a = F{e, D)de/e, 
where 



F{e,D) = 



2[l + (7D + 23)e2] 



7D + 9 
After integration, we obtain 



a — flini exp 



{7D + 23)(e^ - ep„,) + 21og (e/e.ni) 
7D + 9 



(4) 



(5) 
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Figure 1. Schematic illustration of tidal interaction in the system. 



where the subscript ini indicates initial values. It should be empha- 
sized that, when stellar tide is neglected {D — oo) we have 



: aini exp(e 



(6) 



indicating that the final position of the planet is completely deter- 
mined by initial values of the elements. 

It is important to note from equations ([T) and (|2) that the cu- 
mulative effect of tides produces orbital decay and circularization 
on time-scales which varies according to ki /Qi - values (see Jack- 
son et al. 2008; Rodriguez & Ferraz-Mello 2010). Moreover, the 
stellar tide is the only source of orbital decay after total circular- 
ization is achieved. However, we see from equation (|3) that the 
strength of the stellar tide is proportional to the planetary mass mi. 
Hence, tides on the star due to large planets are more efficient than 
tides produced by small planets, in order to induce a planetary mi- 
gration on time-scales comparable to the age of the system. 



3 NUMERICAL SIMULATIONS OF THE TIDAL 
EVOLUTION OF TWO-PLANET SYSTEMS 

In this section we study the dynamical evolution of a system com- 
posed by two planets orbiting a central star. The purpose is to in- 
vestigate how the presence of an outer companion affects the tidal 
evolution of the inner planet. 



3.1 Numerical model 

We will consider systems in which the inner planet is a super-Earth. 
As discussed in the previous section, in this case, the contribution 
due to stellar tides is negligible due to the small planet mass. We 
assume that both planets are deformed under the tides raised by 
the central star. In addition to tidal force, we consider that the in- 
ner planet is also affected by the correction of Newtonian potential 
of the star due to general relativity. The reference frame chosen 
is centered in the star and the motion of the planets occurs in the 
reference plane (i.e coplanar motion). 

According to our assumptions, the orbits evolve under the 
combined effects of mutual interaction and tides raised by the star. 
The equations of motion of the planets are written as 



ri 



G(mo + mi] 



ri + Gm2 



V2 - ri 
ir2 - ri|^ 



r2 



(V) 



momi mo 



V2 = 



G(mo + m2) 



r2 + Gmi 



ri - r2 

|ri — r2p 



(mo + m2) (f + frci) 



mom2 



mo 



(8) 



Table 1. The adopted orbital elements and physical data of the CoRoT-7 
system (Leger et al. 2009; Queloz et al. 2009, Ferraz-Mello et al. 201 1). The 
value of R2 was computed assuming an Earth mean density for CoRoT-7c. 
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where frd is the general relativity contribution acting on the inner 
planet, which is approximately given by 



frel 



Gmimo 



Vi ) ri + 4(ri ■ vi)vi 



(9) 



where vi = ri and c is the speed of light (see Beutler 2005). f and 
g are the tidal forces acting on the masses mi and m2, respectively. 
We use the expression for tidal forces given by Mignard (1979): 



f = -SfciAti- 



-[2ri(ri • vi) +r^(ri X fii +vi)], (10) 



= -3fc2Afr 



[2r2(r2 ■ V2) + r2(r2 X n2 + V2)], (11) 



where fli is the rotation angular velocity of the i-planet, for i — 
1, 2. It is worth noting that the Mignard's force is given by a closed 
formula and, therefore, is valid for any value of eccentricit}0. Ati 
is the time lag and can be interpreted as a delay in the deformation 
of the tidally affected body due to its internal viscosity. 

The time lag is related to the dissipation function (Q) of the 
deformed planet, which is a quantity used in the study of tidal evo- 
lution of extra-solar planets. However, in many classical theories, 
the tidal potential appears written as a function of phase lags (ej) of 
tidal waves, introduced in each periodic term in order to consider 
the effect of viscosity (Darwin 1880; Kaula 1964; Ferraz-Mello et 
al. 2008). In fact, the expansion of equations dlOt-dl lb leads to the 
same forces used in the Darwin's theory when the phase lags of the 
tide components are assumed proportional to their frequencies (the 
so-called linear model): Sj = vjAt, with the same At for all j- 
frequencies. For small lags, Q can be associated to the phase lag of 
the tidal wave through e = 1 /Q (see Efroimsky and Lainey 2007, 
Efroimsky and Williams 2009 for a complete discussion). For plan- 
ets in stationary or synchronous rotation, the tidal evolution is well 
described by the phase lag corresponding to the frequency n (an- 
nual tide, see Ferraz-Mello et al. 2008). Thus, in the linear model 
we have l/Q ~ nAt or, using a modified definition which absorbs 
the Love number, Q' = 3Q/2k, we obtain kAt = 3/{2Q'n) (see 
Leconte et al. 2010). 



3.2 Application #1: CoRoT-7 planets 

We numerically investigate the recently discovered CoRoT-7 plan- 
etary system, which is composed by two short-period planets or- 
biting a central star. The star, CoRoT-7, is a G9V with mass 
mo = 0.93 mQ, radius _Ro = 0.87i?Q and age of 1.2-2.3 Gyr 
(Bruntt et al. 2010). The inner planet, CoRoT-7b, and outer planet. 

Note however that, for a more accurate description, a large number of 
harmonics in the expansion of the tidal potential should be considered when 
the star-planet distance is small enough (see Taylor and Margot 2010) 
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Figure 2. Time variation of semi-major axes for the planets of CoRoT-7 
system. The excitation of ei enables the orbital decay of CoRoT-7b due 
to tides on the planet raised by the star. CoRoT-7c also migrates due to 
tides but with a slower rate, resulting in a divergent migration between the 
planets. 



CoRoT-7c, are assumed to have the masses given in Ferraz-Mello 
et al. (2011), mi = 8.0 and m2 = 13.6 mg, respectively. 
The radius of CoRoT-7c is not known. We adopted values of Q' 
and mean density consistent with a rocky terrestrial planet. Their 
orbital periods are PioA = 0.854 and P2orb = 3.698 days (Leger 
et al. 2009; Queloz et al. 2009). CoRoT-7b is the third exoplanet 
with shortest orbital period discovered up to the present date. 

The adopted physical and orbital parameters of the system are 
listed in Table [T] As we will shown, in the current configuration 
of the system, both orbits must be circular. We emphasize that the 
masses of the CoRoT-7 planets satisfy the relation vnxjra-i < 1. 
Choosing Q'l = Q2 ~ 100 (a typical value for terrestrial planets) 
and using current values of semi-major axes, we obtain kiAti = 3 
min and fc2At2 ~ 12 min. A scaling factor is often used in the 
numerical simulations in order to accelerate the process of tidal 
evolution. In the simulations reported in this paper, we did not use 
a scaling factor. We also assume that the inclinations of both orbital 
planes with respect to the reference plane are zero. 

The next step to perform the simulation is the choice of the ini- 
tial configuration of the CoRoT-7 system. (It should be emphasized 
that, since we are studying a dissipative process, we cannot use a 
simple backward simulation for that sake). The choice is arbitrary, 
but it can be constrained by some considerations. We assume that 
the planets have experienced an orbital decay due to tidal effects, 
thus its past location must be more distant from the central star than 
the current one. Finally, as will be discuss in Section |4] the angular 
momentum of the system can be considered as invariable during 
the evolution of the system under the interplay of the gravitational 
and tidal forces. Gathering this information, we assume the follow- 
ing initial values for the orbital elements: enni — 0, e2im ~ 0.2, 
aiini = 0.0188 au and a2ini = 0.0462 au. 

Figs. |2] and [3] show the result of a numerical integration of 
equations l|7j and l[8j using the RA15 code (Everhart 1985). It is 
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Figure 3. Time variation of eccentricities and Aro (in degrees) for the plan- 
ets of the CoRoT-7 system. The eccentricity of the initial inner circular orbit 
is excited by the outer companion. Tides on the inner planet circularize its 
orbit in approximately 7 Myr. The eccentricity of the outer orbit is also 
damped because of the total angular momentum conservation (see Section 
|4). At the end of the simulation, a state of double circularization is obtained 
(see text for discussion). 



shown the orbital decay and circularization of the orbits due to 
tidal effects. The inner planet eccentricity initially undergoes large 
variations, but, after a short time of damped oscillations, it reaches 
a quasi-equilibrium value (Mardling 2007, see Sec. 13.2. It . After 
relaxation, ei decreases smoothly to reach near-circularization in 
about 7 Myr. The eccentricity of the outer planet is also damped 
and a pair of circular orbits is obtained as final configuration of the 
system. 

The excitation of ei activates the orbital decay of the inner 
planet, which requires a non circular orbit according to equation 
([TJ (with s — 0, because we are not considering the contribution of 
stellar tides in our simulations). The migration of CoRoT-7b toward 
the star occurs until the circularization of its orbit. Meanwhile, 02 
also decreases (and is only weakly affected when the 4/1 mean- 
motion resonance is crossed). However, when equations (llO) -lll It 
are expanded up to third order in eccentricities (see Ferraz-Mello et 
al. 2008), we obtain that that the ratio of forces on CoRoT-7b and 
CoRoT-7c is approximately 16, indicating a more intense migration 
of the inner planet. 

The bottom panel of Fig.[3]shows the time variation of the sec- 
ular angle Aw = wi — zu2, where zu means longitude of pericen- 
ter. Initially, Azu oscillates around Aw — 0° with large amplitude; 
this stage of the evolution corresponds to the damped oscillations 
of the planetary eccentricities. Reaching the relaxed state, in which 
the system oscillates with amplitudes close to zero, the oscillation 
amplitude of the secular angle also tends to zero. 

It is worth mentioning that, due to the divergent migration be- 
tween the planets, the capture in mean motion resonance is not pos- 
sible. 
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Figure 4. Quasi-equilibrium e^'' given by equation {12), as a function of Figure 5. Time variation of e2 in the CoRoT-7 system. Tlie comparison is 

the mass ratio. The vertical dashed line corresponds to the mass ratio of done between the integration of equation (B] and the result of numerical 

the CoRoT-7 system. The curves show the cases with and without general simulation (see text for more discussion), 
relativity. 



3.2.1 The short-time evolution of the eccentricity of CoRoT-7h 

Fig. [3] shows that, approximately at t — 0.5 Myr, the amplitudes 
of oscillation of the eccentricities are damped and the variation be- 
comes smooth. This situation is called quasi-equilibrium. Mardling 
(2007) derived an expression to compute the quasi-equilibrium 
value of the inner planet eccentricity, e^'', using the Legendre ex- 
pansion of the disturbing function, up to fourth order in ai/a2. It 
is given by 



(5/4)(ai/a2)e2e2 



1 - \/a^Ja^{mi/rn.2)£2 ^ + Tell 



(12) 



where £2 = yl — ef and 7 = 4(niai/c)^(mo/m2)(a2/ai)^. 

Fig. |4] shows the variation of e^'' given by equation J12b . for 
CoRoT-7, as a function of mi/m2. The relevant value of 62 which 
should be used is the one when e^'' is reached. In addition, the con- 
dition l ll2b was obtained under assumption that semi-major axes are 
constants. We show the cases with (7 7^ 0) and without (7 = 0) the 
contribution of general relativity. The inclusion of 7 7^ reduces 
the value of e^''. Hence, since the mean rate of orbital decay is pro- 
portional to e\, the omission of the relativistic potential can result 
in an underestimation of the circularization and orbital decay time- 
scales. Mardling & Lin (2004) also showed that the inclusion of the 
relativistic potential is very important when the planet approaches 
the star. 



62 



A 



£2 



(13) 



where rdrc is the circularization time-scale of the inner orbit in the 
single-planet case (i.e. ei/ei, where e\ is given by equation l[2}); 
A = (25/16)(mi/m2)(ai/a2)^/^ and F{e2) = el{l - aej ' + 
762)^, with a = (mi/m2)(ai/a2)^''^. 

The comparison between the integration of the equation ( I13t 
and the result of the simulation of exact equations of motion {T} 
and l|8) is shown in Fig. |5] The difference between both results 
may arise from the fact that for obtaining equation (113) is explicitly 
assumed that semi-major axes are constant (see Mardling 2007]j. 
In Section [33] we will show that the disagreement increases in a 
system for which m\/m2 1. 

We note that, according to equation J13t . the rate of €2 - 
damping decreases with 7712 /mi, enabling a non-circular outer or- 
bit be maintained over long times. In addition, e^"* increases with 
62 as equation l ll2t indicates. Hence, since the mean rate of orbital 
decay is proportional to ef, a significant inner orbital shrinkage can 
occur for large m2/m\ (see next section). 



3.2.3 The rotation of CoRoT-7b and CoRoT-7c 

The tidal force produces a torque that directly affects the rotation 
evolution of the deformed bodies, in such a way that Cif2i = ri x 
f and C2ft2 = r2 X g. Thus, 



+ r, X Vi], 



(14) 



3.2.2 The long-time evolution of the eccentricity of CoRoT-7c 

Our numerical simulation has shown that the eccentricity of the 
outer planet is strongly affected during the orbital decay (see Fig. 
[2I top panel). However, the source of 62 - damping is not the direct 
action of tides on the outer planet since, as we have already shown, 
the magnitude of the tidal force on CoRoT-7c is much smaller than 
the corresponding force on CoRoT-7b[f| 

Mardling (2007) also studied the long-term variation of the 
outer planet eccentricity and provided an expression to calculate 
the 62 - damping: 



where d is the polar moment of inertia of the i-planet. 

To obtain the time evolution of the angular velocity Q,i, the 
above equations (i = 1,2) were integrated simultaneously with 
the equations of motion l[7j-l[8]l. The result of the integration is 
shown in Fig.|6] The rotation periods, defined as Pi rot = 2% /Q,i, 
reach stationary values, which, up to second order in 6^, are given 
by Pirot = Piorb/(l + 6ef) (Hut 1981; also see Ferraz-Mello et 
al. 2008; Correia, Levrard & Laskar 2008). In the same figure, we 
also plot the orbital periods defined as Pi orb = 2n /ui, which de- 
crease as a consequence of the orbital decay. The synchronization 



* One can take an idea about which is the contribution of tides raised in 
the outer planet to the circularization process of its orbit. For that sake, we 
note that the ratio between the circulaiization time-scales associated with 
the direct tidal interaction with the star and the one given by equation (T3) 
is of the order of (02 /ai )^ ^ 1- 



^ We stress that equation l |13) gives an estimate for the damping time-scale 
of 62 (the full expression is given in equation (57) of Mardling 2007). How- 
ever, it is shown in that work that the approximate solution ( |13) underesti- 
mates the time-scale on which 62 evolves. Hence, the difference between 
analytical and numerical results would be even large if the full expression 
were used. 
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Figure 6. Time variation of the rotation and orbital periods of CoRoT-7 
planets. Synchronization is reached as the doubly circular state approaches. 



0.2 
0.15 

0.1 
0.05 



— 1 — I I 1 1 iii| 1 — n-mTTjT — I — I I 1 1 1 1 

y=0 





0.01 




_l ' 



0.1 



10 



m.|/m2 



Figure 7. Quasi-equilibrium e^'' in the second application, given by equa- 
tion I12t . as a function of the mass ratio. The contribution of general rela- 
tivity does not affect the value of and thus can be neglected. 



Table 2. Orbital elements and physical data coiTesponding to an hypotheti- 
cal Sun-super-Earth- Jupiter system. The value of R\ was computed assum- 
ing a ten'estrial mean density. 
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between the orbital motions and the rotation of the planets occurs 
when Q,i = rii or Pii-ot — PioA- In Fig. [6] we can observe that 
this condition is only reached after the circularization of the planet 
orbits. 

3.3 Application #2: super-Earth with a giant companion 

In this section, we investigate the tidal evolution of a system in 
which the inner planet is much smaller than the outer one. We chose 
the fictitious system composed by a super-Earth and a Jupiter-like 
outer companion, with masses mi = Sm® and m2 = Imj, 
respectively (mj is the mass of Jupiter). Note that m\/m2 — 
0.0157 <C 1. The initial configuration and physical parameters of 
this hypothetical system are listed in Table |2] in which the values 
of semi-major axes correspond to orbital periods of 2.92 and 11.6 
days. The large orbital period of the outer planet, in addition to the 
large Q-value for Jupiter-like planets, enable us to neglect the tides 
on the outer companion. Moreover, as shown in Fig.|7]for the con- 
sidered system, the contribution of the relativistic potential is also 
negligible. Here, we also assume zero inclinations. In addition, the 
calculation of fciAti = 3/(2Q'ini) with Q'^ = 100 results in 
fci Ail = 10 min. 




20 30 40 
time [Myr] 



0.12 




20 30 
time [Myr] 

Figure 8. Time variations of orbital elements in the case mi/m2 <C 1 
(c.f Table|2). The evolution in the case in which the super-Earth is the only 
planet is also shown. The inner orbit reaches a quasi equilibrium value in 25 
Myr. retarding the orbital circularization and enhancing the orbital decay. 



outer companion delays the circularization of the inner planet orbit 
(see Fig.|9]l. As a consequence, the presence of the second planet 
lead to a faster decay of the inner planet orbit. (In the single-planet 
case, the semi-major axis a\ decreases only slightly and reaches a 
final constant value (0.0396 au) after total circularization^l 



3.3.1 Short-time evolution (50 Myr) 

Fig.[8]shows the time evolution of the inner planet semi-major axis 
(solid curve in the top panel) and two eccentricities (bottom panel) 
over the first 50 Myr. Due to the large mass of the outer planet, 
62 is almost constant over the time interval shown in the figure (see 
equation l|13ll). On the contrary, ei initially suffers large oscillations 
which are quickly damped to the quasi-equilibrium value. On the 
same graph, we plot the variation of ei in the case of a single- 
planet system, when the interactions with one outer planet do not 
exist. Comparing the two cases, we note that the presence of the 



3.3.2 Long-time evolution (300 Myr) 

Fig.|9]shows the time variation of the planet semi-major axes (top 
panel) and the eccentricities (bottom panel) over 300 Myr. After 
e^'' is reached, the inner planet eccentricity decreases smoothly to 
very small values. The outer eccentricity follows the evolution of 
the inner one and, at the end of simulation, about ~ 0.25 Gyr, a 
pair of circular orbits is obtained. When the time approximates to 

^ This value can be obtained using equation ^6) with e = 0, ejni = 0.1 and 
ami = 0.04 au. 
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Figure 9. Long-term evolution of semi-major axes and eccentricities in a 
Sun-super-Earth- Jupiter system. Due to orbital decay, the final position of 
the inner planet is close to the Roche limit. The value of ei is damped to 
zero and the outer orbit is also circularized. 



0.12 Gyr, ei is strongly excited due to the passage of the system 
through the 5/1 mean-motion resonance but, immediately after that, 
is damped to its previous value. 

The orbital decay of the inner planet occurs until the to- 
tal circularization of its orbit (see top panel in Fig. |9). Note that 
the final position of the super-Earth is close to the Roche limit 
fflRoche = (7?i/0.462)(mo/mi)^''^ (Faber, Rasio & Willems 2005), 
which gives ORochc = 0.0064 au. The outer planet semi-major axis 
does not varies in the considered time interval. 

In Fig. [To] we compare the damping of 62 resulting from nu- 
merical simulation (dashed line) and that obtained through the in- 
tegration of equation ( |13t of Mardling's model (solid line). The re- 
sults are very different, indicating that equation (1131 is not a good 
approximation to compute the decreasing of 62 in the case of the 
massive outer companion. The discrepancy can be understood tak- 
ing into account that, for this particular system, there is a signifi- 
cant orbital decay due to tides on the inner planet, while the equa- 
tion ( |13t was obtained assuming A and a as constants, i.e. constant 
semi-major axes (Mardling 2007). Thus, it is shown that the process 
of outer planet circularization is accelerated in the case of signifi- 
cant orbital decay of the inner planet, if compared with the result 
of equation il3t . 



4 CONSERVATION OF THE ANGULAR MOMENTUM 

Many dynamical features of the planet evolution shown in the pre- 
vious sections, can be explained considering the conservation of the 
angular momentum of the system. The total angular momentum of 
the three-body system in the astrocentric reference is the vectorial 
sum of the orbital momenta of the planets, L°* (i = 1, 2), and the 
rotational momenta of the planets and star, L'^"' (i = 0, 1, 2). 

If Vi is the astrocentrical position of the i-planet, the orbital 
component of the angular momentum is given by 




0.1 1 10 

Time [Gyr] 

Figure 10. Time vaiiation of 62 in a Sun-super-Eaith- Jupiter system. The 
comparison is done between the integration of equation 03) and the result 
of numerical simulation (see text for more discussion). 



Lorb — rrnri ^ ~ rmri x rmi- 



(15) 



where AI = 5^^_q rrii. Note that the second term of the equation 
l |15l l appears due to the use of the astrocentric reference frame. This 
term is of second order in masses and, assuming rrii <^ mo, can 
be neglected. Moreover, x = G(mo + mi)ai{l — el) k, 
where k is a unitary vector normal to orbital planes. Hence 



ik, 



(16) 



where m[ = mi^Gmo. 

The rotation component of the angular momentum is given by 



E 

i=0 



where d are the moment of inertia with respect to the axes of ro- 
tation and fli are the angular velocities of rotation, for i = 0, 1, 2. 
Because we account only for the planetary tides on the planets, fli 
and fl2 vary with time, while flo is constant. 

We can show that for planets, L"' < Lf*"; thus, the contri- 
bution of the rotational component to the total angular momen- 
tum can be neglected. Indeed, using the third Kepler's law, we 
have Lf /Lf = ^,{R^/a,f{n^/n^), where we have introduced 
Ci = ^iiriiR^, with < .fi ^ 2/5. In the case of the quasi- 
synchronous planet discussed in Section [3.2.3l when Q.i ~ n^, this 
ratio is small, of the order of {Ri/ui)'^. 

In contrast, for the star, the quantity / Lf^ = 
io{'mo/mi){Ro/aif'{yio/ni) is close to unity, when we consider 
close-in planets and a slow-rotating star, for which the condition 
fio ^ rii is satisfied. Hence, the total angular momentum of the 
system can be approximately written as 



L ~ {m\^ai{l - ei) + mj V«2(l - e^) + CoQo)\i. (17) 

Is should be noted that some quantity of angular momentum 
is removed due to the orbital decay of the inner planet. Thus, ac- 
cording to above equation, for small ei, the eccentricity of the outer 
planet would be damped in order to preserve the angular momen- 
tum conservatation. 



4.1 The orbital decay 

As a consequence of the total angular momentum conservation, we 
can obtain the minimal possible value for oi / 02 during the orbital 



8 A. Rodriguez et al. 



0.04 



S- 0.03 



0.02 



0.01 







- eiini=0.1 




e2ini=0.1 




- aij„j=0.04au 




a2ini=0.10 au 





10 
m2/m^ 



100 




Figure 11. Plot of equation )19t showing the dependence of the minimum 
value of the inner planet semi-major axis with m2 /rrti . The initial elements 
of Table |2] were used to compute ki and K2. The case of a more massive 
inner planet is shown in the bottom panel. 



decay, knowing that this minimum value should be achieved once 
orbital circularization is reached. Imposing ei = 62 = in equa- 
tion l |17t we obtain 



aiir 



02 



= (i' - 1)^ 



/m2 
Vmi 



(18) 



where L' = {L — Co^o) /m^y/oa and m'2/m'i ~ m2/m\. For 
those cases in which tides on the outer planet can be neglected, 02 
and thus L' are constants of the problem. Applying equation ( 118) 
to the super-Earth- Jupiter system we obtain ai„,in/a2 — 0.096, 
reproducing the corresponding result of the numerical simulation. 
Hence, it is worth noting that the final position of the inner planet 
can be obtained from the analysis of the angular momentum con- 
servation, as given by equation dlSt . 

Equation l |18b allows us to investigate the variation of aimin/a2 
for several mass ratios. However, we note that there is also an im- 
plicit dependence on m2/mi through L', which is not shown in 
equation dlSt . Hence, in order to express the full dependence we 
define two constants, namely, k\ = (aimi/a2ini)(l — efini) and 
'^2 = A^/l — e2iiii- Using the definition of L' and replacing it into 
equation l llSt , we obtain 



(K2 — 1) h Kl 

mi 



(19) 



Fig. [TT] shows the variation of aimi„ with the mass ratio 
1712 /m-i according to equation ( 119) , where we have used initial val- 
ues of the orbital elements listed in Table |2] We see that, for fixed 
initial values of the elements, the inner planet will come closer to 
the star for high-mass companions. Therefore, as discussed in Sec. 
13.2.21 the orbital decay of the inner planet is enhanced. 

For sake of completeness, we show in the bottom panel of Fig. 
II lithe result for 1x12 /m\ < 1, which is the case of a more massive 
inner planet. Note that the orbital decay is very weak, supporting 
the assumptions done in Mardling (2007) of constant semi-major 
axes in the case of slow-mass outer companions. 



5 EVOLUTIONARY ROUTES IN THE PHASE SPACE 

In this section we will show that evolutionary routes of the migrat- 
ing system are traced by stationary solutions of the conservative 
secular problem. For this task, we explore the secular dynamics of 
the planar two-planet system, described in detail in Michtchenko & 
Ferraz-Mello (2001) and Michtchenko & Malhotra (2004). 

In the context of hamiltonian formalism, the stationary val- 
ues are defined by the secular Hamiltonian function, in which the 
fast-period terms, involving arguments with combinations of mean 
longitudes, are removed by first-order averaging. In this case, the 
secular problem of coplanar two-planet systems is reduced to one 
degree of freedom. Regarding to the angle Avj, the stationary so- 
lutions are restricted to Azu = 0° and Avj = 180°, which corre- 
spond to Mode I and Mode II solutions, respectively. In addition, 
the stationary values of eccentricities for these modes are e*j and 
e*ji, for i = 1, 2. Periodic solutions of the secular system are cir- 
culations or oscillations around Modes I or II and the stationary 
solution can be represented by a point in the plane of eccentricities. 
The position of this pont in the phase space depends on the pa- 
rameters of the secular problem, namely, angular momentum and 
semi-major axes ratio. 

To understand how the stationary solutions are affected by the 
introduction of dissipative forces on the system, we will assume 
that the dissipation is very slow. This means that the dissipation rate 
is much smaller than the characteristic time of the secular system, 
which is the proper period of the secular angle Azu. Under this 
assumption, the motion of the pair of planets is nearly conservative 
and, for very short time intervals it can be described by the classical 
secular theory. 

For the sake of simplicity, let us consider that periodic solu- 
tions are oscillations around Mode I. When the tidal effects are 
taken into account, the semi-major axis of the inner planet de- 
creases slowly. Since the semi-major axes ratio is changed, the po- 
sition of the equilibrium (e^j, e^/) in the phase space will be also 
changed, in order to preserve the total angular momentum of the 
system. In such a way, the system oscillates around a centre whose 
position in the phase space changes slowly due to dissipation. 

Assuming continuous values of the inner semi-major axis, we 
obtain the equilibria (e^j,e2/), for all possible values of ai. A 
curve in the space (ei, 62), which is composed by successive equi- 
librium points, is a locus of stationary solutions in the space of 
eccentricities and we refer to it as LSE. By definition, all solution 
in LSE curve have the same value of angular momentum. 

Fig.[T2]illustrates the LSE curve for Mode I in the case of the 
CoRoT-7 planetfl comparing with the curve resulting from numer- 
ical simulation of the system. For sake of simplicity, we do not in- 
clude the relativistic potential into the analysis of mean equations. 
The initial conditions for the numerical integration were chosen as 
Azu = 0, ei = and 62 — 0.2. Initially, the eccentricities os- 
cillate rapidly with decreasing amplitudes around an equilibrium 
which slides slowly along the LSE curve. When ei ~ 0.05 and 
62 — 0.08, the amplitude of oscillation becomes close to zero and 
the solution is almost coincident with the LSE curve, until the sys- 
tem reaches the total orbital circularization. 

Detailed discussion about the secular dynamics of non- 
conservative two-planet systems can be found in Michtchenko & 
Rodriguez (2011). 



^ We have seen in Section [X2l that the angle Azu oscillates with almost 
zero amplitude around Aro = 0°, indicating motion around Mode I. 
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Figure 12. Compaiison between computation of LSE curve and the result 
of numerical simulation for CoRoT-7 system (without the contiibution of 
the relativistic potential). After the angle Aro is close to 0° (Mode I), both 
curves are almost coincident. 



Figure 13. The mean hfetime Ta ^ of CoRoT-7 planets due to stellar tides 
in function of the dissipation on the star. For ko/Qo = 10^'', Ta^ = 0.97 
Gyr and = 367 Gyr 



6 THE ROLE OF STELLAR TIDES 

In this section, we evaluate the time required for CoRoT-7 planets 
reach the Roche limit with the star. As seen in Section|2] in the case 
of small-mass planets, the contribution of tides raised by the planet 
on the star can be neglected if compared with the contribution of 
tides raised by the star on the planet. However, we know that after 
orbital circularization, the stellar tide is the only source of orbital 
decay. Thus, due to the current orbital configuration of the CoRoT-7 
system, with two circular orbits, only the stellar tides could produce 
the orbital decay. 

We define the mean lifetime of CoRoT-7b as the time for 
which the Roche limit is reached, starting with ai = aicuncnt = 
0.017 au. According to equation ([T) with p = 0, ei = and in the 
limit mi ^ mo, the mean rate of semi-major axis variation is 



3\/Gmo ko nil 5 
Equation i20\ can be integrated as 

r "Roche 

J 11/2 

dai — —ItTai, 



(20) 



(21) 



^ 1 current 

where JC = 3y/Gmo{ko/Qo){'mi/mo)RQ, whereas Ta^ is the 
mean lifetime of CoRoT-7b and aRoche ~ 0.0049 au (~ 1.2 i?o). 
Solving for Tai , we obtain 



„ / 13/2 13/2 ^ „ 

2 I'^l current ~ %ochc J ^0 "^0 



(22) 



39 ^/GW7^g ko mi ' 

Fig. [H] shows the variation of r^^ as function of ko/Qo ac- 
cording to equation ( 122) . It is easy to see that the mean lifetime of 
CoRoT-7b is small for high values of ko/Qo- On the other hand, 
increases when the energy within the star is not efficiently dis- 
sipated (small ko/Qo )■ For a typical value of ko/Qo = 10~®, 
we obtain Ta^ = 0.97 Gyr (see also Jackson, Barnes & Green- 
berg 2009). If the same analysis is done for CoRoT-7c, we obtain 
flRoche = 0.0062 au (~ 1.5 Ro) and Ta-^ = 367 Gyr, indicating that 
CoRoT-7c should not be destroyed by tides within the lifetime of 
the star. 



7 DISCUSSION AND CONCLUSIONS 

In this work we have investigated the tidal evolution of super-Earth 
planets with exterior companions. Previous works have already dis- 
cussed that problem, giving special attention to the case of a more 



massive inner planet. In Mardling (2007), explicit equations were 
obtained giving the variation of eccentricities in the context of con- 
stant semi-major axes (see equations J12t and (I13t). Here, we cen- 
tre the attention to the case of a more massive outer companion. 

Through an approach in which the exact equations of motions 
are numerically solved, we have shown that the super-Earth orbital 
decay is accompanied by orbital circularization of both planets. In 
addition, we show that the rate of e2-damping is small for large 
values of m2/mi, enabling to sustain a non-circular outer orbit 
for longer times. Thus, as e'^ increases with 62, a significant inner 
orbital decay can occur for large ra2/m\, since its mean rate is 
proportional to e\. 

As a real application of tidal evolution of a super-Earth accom- 
panied by a more massive outer planet, we investigated the CoRoT- 
7 two-planet system. We have shown that, starting from an initial 
circular orbit, ei is excited enabling the activation of the orbital 
decay of CoRoT-7b (which requires a non-zero value of ei). The 
equilibrium value of ei is quickly damped by tides on the planet 
and a strong decreasing of a\ is observed in the short-term tidal 
evolution. 

In addition, the long-term tidal evolution has shown that the 
current orbital configuration of the system is achieved, including 
the double-circularization state of CoRoT-7 planets. After circular- 
ization, the migration toward the star continues due to the action of 
stellar tides. We have shown that CoRoT-7b can reach the Roche 
limit in approximately 1 Gyr for k\/Q\ = 10~^ (see also Jackson 
et al. 2009). 

The study of stationary solutions of mean equations is useful 
to better understand the secular dynamics in the space cos I\vj - 
ei sin The periodic solutions around Mode I (aligned pericen- 
ters) and Mode II (anti-aligned pericenters) may be easily identified 
in the conservative case. In addition, the analysis can be extended 
to the case in which dissipation is included. It is possible to con- 
struct a curve in the space of eccentricities composed by successive 
stationary solutions, which are obtained one by one through the 
same technique used in the conservative case. We have shown that 
the result of numerical simulation of exact equations for CoRoT-7 
planets is in good agreement with the curve of stationary solutions 
of mean equations of motion. 

Future studies involving the dynamical evolution of a two- 
planet pair including tidal dissipation will help to investigate the 
past history of the CoRoT-7 system (i.e formation and subsequent 
evolution). This may bring more information about how the current 
orbital configuration of the system was achieved. 

As an example in which m\/m2 <^ 1, we have also investi- 
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gated the tidal evolution of a Sun-super-Earth-Jupiter-like system. 
The results of numerical simulation have shown that 62 is not well 
described by integration of the simple model represented by equa- 
tion J13b . which describes the time variation of 62 in the long-term 
evolution. We have shown that the doubly circular state is acceler- 
ated when there is a significant orbital migration of the inner planet, 
showing that the final configuration of the system is affected by the 
shrink of the inner orbit. 

A simple analysis based on the angular momentum conserva- 
tion shows that it is possible to predict the final value of the inner 
planet semi-major axis, provided the semi-major axis of the outer 
planet is fixed. 
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